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Abstract 



Q^ ' The angular intensity distribution of He beams scattered from compact clus- 

^ I ters and from diffusion limited aggregates, epitaxially grown on metal sur- 



faces, is investigated theoretically. The purpose is twofold: to distinguish 
compact cluster structures from diffusion limited aggregates, and to find ob- 
servable signatures that can characterize the compact clusters at the atomic 
' level of detail. To simplify the collision dynamics, the study is carried out 

in the framework of the sudden approximation, which assumes that momen- 
tum changes perpendicular to the surface are large compared with momentum 
transfer due to surface corrugation. The diffusion limited aggregates on which 
the scattering calculations were done, were generated by kinetic Monte Carlo 
simulations. It is demonstrated, by focusing on the example of compact Pt 
Heptamers, that signatures of structure of compact clusters may indeed be 
extracted from the scattering distribution. These signatures enable both an 
experimental distinction between diffusion limited aggregates and compact 
clusters, and a determination of the cluster structure. The characteristics 
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comprising the signatures are, to varying degrees, the Rainbow, Praunhofer, 
specular and constructive interference peaks, all seen in the intensity distri- 
bution. It is also shown, how the distribution of adsorbate heights above the 
metal surface can be obtained by an analysis of the specular peak attenuation. 
The results contribute to establishing He scattering as a powerful tool in the 
investigation of surface disorder and epitaxial growth on surfaces, alongside 
with STM. 
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I. INTRODUCTION 



The process of epitaxial growth of metal or semiconductor films on a surface presents 
some major theoretical, experimental and technological challenges. In epitaxial growth, 
metal or semiconductor atoms are adsorbed on a corresponding surface under thermal condi- 
tions, to form two- and three-dimensional structures on top of it. The physical and chemical 
properties are determined by the final form of these structures. These may be of dramatic 
importance, e.g, in the production of electronic devices. One of the most exciting aspects 
of epitaxial growth kinetics, is that it prepares disordered structures in the interme- 
diate stages. The disorder manifests itself in the formation of various types of clusters or 
diffusion-limited-aggregates (DLA) on top of the surface. These structures may be monolay- 
ers (usually at high temperatures, when the diffusivity is large), in which case the disorder 
is two-dimensional, or they may be composed of several layers, giving rise to disorder in 
three dimensions. Epitaxially grown structures of this type have for some time now received 
wide attention, as they offer an exceptional opportunity for both experimental and theoret- 
ical study of disorder. Clusters in particular play an important role in the epitaxial-growth 
kinetics, as it is the continual process of their formation, restructuring and dissociation, 
through which growth actually takes place. Issues such as the structure, stability and sur- 
face diffusion of metal clusters have been addressed experimentally by Wang and EhrlichS. 
Rosenfeld et al.i have recently suggested the intriguing possibility of the formation over 
a certain temperature and coverage range, of a high density of small, stable and compact 
clusters of seven ad-atoms of Pt (Heptamers), during the process of epitaxial growth of Pt 
on Pt(lll). A recent theoretical first-principles study by Stumpf and Scheffleri predicts 
hexagonally shaped islands for Al on Al(lll). These developments tentatively suggest that 
metastable states of systems consisting of clusters of well-defined shape, and perhaps size 
may exist. Regarding DLA, several authorsi^ have presented simulations which reveal the 
complexity anticipated in the shapes of these structures, in which no regular clusters can 
be discerned. As these and other studies show, a large variety of structures contribute to 
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surface disorder. No satisfactory and comprehensive theory of the epitaxial growth process 
is as of yet available, much due to the absence of reliable interaction potentials between 
the ad-atoms and the surface. Progress at this stage thus hinges critically on experiments. 
To a large extent the information obtained from different modern methods is complemen- 
tary. Since the various techniques probe different aspects and distance scales of the surface 
structures, the most comprehensive picture will probably be attained via an integration of 
results from all techniques. TEAS and in particular He-beams, has already proved among 
the most useful tools in analyzing surface phenomena. Experiments by Poelsema and co- 
workersiH0 and Ernst and coworkersliUli have employed He-scattering in dealing with ques- 
tions concerned with scattering from isolated surface defects, and were so able to extract 
a wealth of information on the defects and their mutual attractions. The attenuation of 
specular intensities, in particular, proves an extremely useful tool in these investigations. 
Theoretical studies by Gerber and co-workersEUll utilizing both time independent (Sudden 
approximation) and time-dependent methods (quantum wave-packets), have analyzed the 
calculated angular distribution and specular peak attenuation. These authors have been able 
to identify connections between the features of the distribution and the attenuation, and 
microscopic characteristics related to the surface defects. An important advantage offered 
by He-scattering, is that He is capable of probing structure both at the local and global 
levels. The interpretation of He-scattering experiments, however, is rather more involved 
than that of the direct-imaging techniques such as FIM and STM. The purpose of the present 
article is to address the questions of whether "signatures" of the various shapes assumed by 
the clusters and DLA are present in the spectra obtained with He-scattering experiments, 
and to find how they may be extracted with the proper analysis. The answers to these 
questions are crucial to the understanding of epitaxial-growth induced surface disorder, and 
to enhance the usefulness of TEAS as an indispensable tool in the study of such surface 
phenomena. The experiments by Rosenfeld et ahS and Bott et al.i^l on Ft and by Ehrlich 
et al.0 on Ir suggest that under proper conditions compact clusters may be observed as a 
dominant phase, whereas DLA will predominate at other conditionJ^^. This is supported by 



Stumpf and Scheffler's theoretical calculations for Aid. It will be our main purpose in this 
article to demonstrate that compact clusters can easily be distinguished from DLA in a He- 
scattering experiment. Furthermore we hope to demonstrate, using the example of compact 
Heptamers study for clusters, that specific structure-related characteristics of com- 

pact clusters can indeed be extracted from an analysis of the angular intensity distribution, 
and that these characteristics, when combined, provide an almost unique signature of the 
structure. Thus by a detailed examination of the angular intensity distribution calculated 
from several model systems, we seek to show that surface-adsorbed structures exhibiting var- 
ious types of disorder, can be characterized by means of He-scattering experiments. The 
structure of the article is as follows: In Section |I| we discuss our model systems. Section |H 



briefly discusses the methods employed for studying the scattering intensities. Section |V 
describes the results of scattering calculations from our model systems. Concluding remarks 
are brought in section 0. 



II. THE MODEL SYSTEMS 

The present study deals with the scattering of a He beam from various types of disordered 
structures, considered pertinent to epitaxial growth modes. The possible formation of such 
structures in kinetic simulations of growth is also treated briefly. We shall treat two types 
of structures: 

1. Compact Clusters on a Surface. 

Within this family, disorder will be represented in an increasing degree of complexity, 
by considering three related models: 

(a) A single Heptamer on an otherwise completely flat Pt(lll) surface. 

(b) A collection of Heptamers on the same Pt surface, adsorbed at random positions, 
but having orientations compatible with the underlying hexagonal symmetry. 
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(c) A collection of Heptamers on the same Pt surface, adsorbed both at random 
positions and with random orientations in the surface plane. 

2. Diffusion Limited Aggregates on a Surface. 

Here we considered non-compact structures, grown under kinetic conditions. 

The clusters and the underlying surface were treated as a static, nonvibrating target 
in the scattering calculations. Clearly, vibrations of the adsorbate atoms and the surface 
are expected to affect quantitatively the scattering distributions measured in TEAS experi- 
ments. However, results of calculations using a rigid, nonvibrating surface system should be 
useful at least for studying the main qualitative effects. Moreover, for diffraction scattering 
from crystalline surfaces, the effect of surface vibrations on the scattering intensities can 
be represented approximately by a simple Debye- Waller factoiH. A similar description is 
expected to be successful also for the scattering intensities from clusters on surfaces. In the 
context of studies of the specular attenuation, work by Poelsema, Comsa and coworkers^! has 
demonstrated that treating the surface as static is a very useful approximation in analyzing 
specular scattering from disordered surfaces. The Pt surface in the scattering calculations 
reported here was treated as completely flat. For a very low corrugation surface such as 
Pt(lll), this should be a good approximation. We now describe in some detail the three 
models of compact Heptamer clusters and the model of kinetic growth of DL A on the surface, 
all of which were used for the scattering calculations. 

A. Single Heptamer on a Flat Pt Surface 

In this section we describe a model for a single, symmetric and static Heptamer on a 
smooth metal surface. The parameters were chosen to represent a Pt Heptamer on a Pt (111) 
surface. Although we do not expect a TEAS experiment to ever probe a surface with just a 
single Heptamer adsorbed upon it, the study of the angular intensity distribution produced 
by such an arrangement is crucial in understanding the specific Heptamer related features. 
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As a motivation for our choice of potential, consider the role of the conduction electrons. 
Clearly the Heptamer is immersed in the sea of surface conduction electrons. However, we 
expect the surface to be much smoother than the Heptamer's top, due to a flow of electrons 
from the top to the surface. The extent to which the Heptamer's nuclei lie bare determines 
both the Heptamer's height above the surface and its top's corrugation. These are very 
important parameters in determining the scattering features. A simple, albeit inexact, class 
of potentials which satisfy the requirement of a corrugated Heptamer-top, are those which 
are pairwise additive in the top. We assume the Heptamer to be composed of an isolated 
cluster of seven Pt atoms adsorbed on a flat Pt surface, with six Pt atoms forming a perfect 
hexagon and the seventh atom located at its center. The interaction between the gas phase 
He and the Pt surface with its adsorbed cluster is assumed additive and is represented by 
the following potential function: 
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Viv) = V,{z) + VHe-Pti\r - Til) (1) 
i=l 

where r = {x,y,z), Vsiz) represents the interaction between He and a flat Pt(lll) 
surface, and V/fe_p((|r — ri|) represents the interaction between He and the ith Pt atom in 
the cluster, located at ri. The coordinate z measures the distance of the He atom from the 
surface plane, while (x, y) are the coordinates in parallel to the surface plane. The functional 
form of Vs(2;) was chosen as: 

-D.{^^- ii^) (2) 

which is a surface scattering version of a Lennard- Jones 6-12 potential. The value Ds = 
A.OmeV{= 8.464 ■ 10~^a.u.), given by Harris et ahil was employed, and Zm was taken as 
7. 87 a.u., which gives a steepness typical of atom- atom interactions with metal surfaces. For 
Vne-pt a Lennard- Jones 6 — 12 potential was used: 

^---(0 = ^ - ^ (3) 

We assumed that the cq interaction parameter between He and a heavy atom X scales 
linearly with the polarizability of X. The polarizability, in turn, was estimated from the 



available atomic volume. In our estimate, the cq parameter so obtained is sufficiently reliable 
for the qualitative purposes of this paper. Thus: 



where 5* is the unit-cell area. This equation is based on taking the continuum limit of a 
sum of pairwise He/Pt interactions. It gives a value of C12 = 3.251 ■ 10''a.u. We do not have 
ab-initio calculations available to provide reliable estimates of our additive-potentials model, 
and thus we do not expect the interaction potential eq.(|I]) to be quantitatively correct, but 
we expect this potential to be realistic on a semiquantitative footing. This should suffice for 
the purpose of searching for qualitative effects in the angular distribution. 



The situation we expect to find in a real experiment is a distribution of clusters of 
various sizes and shapes: isolated Heptamers are probably not realized at the coverages of 
interest for epitaxial growth. Consider as a conjecture the existence of Heptamer cluster 
phases. If confirmed after further study, this distribution will exhibit a pronounced peak for 
symmetrical Heptamers. Hence a more realistic treatment of the scattering off Heptamers 
must involve a distribution of Heptamers on a surface. The Heptameric-phase conjecture 
does, however, leave us with a freedom to consider two models of a distribution of adsorbed 
Heptamers. These correspond to different degrees of disorder, and follow from the possibility 
of either a high or low barrier for rotation of a Heptamer-cluster as a whole, as explained 
below. In this section we investigate the scattering from the two corresponding distributions 
of Heptamers, while ignoring the possible existence of other cluster types. We expect the 
other cluster types not to alter the central features of the intensity spectrum, provided 
the distribution is peaked sharply enough about the Heptamers. In this sense the present 



C6(Pt-He) ~ C6(Xe-He) Vol(Pt)/Vol(Xe) = 6.322 a.u. 



The C12 parameter was extracted from eq.(^ by the formula: 





B. Randomly Adsorbed Heptamers 
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framework is an adequate tool for our purpose of providing signatures of Heptamers. In 
the two models we investigated in this section we considered the case of randomly adsorbed 
Heptamers. In both, the center of each Heptamer was chosen randomly from a uniform 
distribution, and then placed on the nearest lattice site. We excluded those Heptamers that 
overlapped with or bordered on other Heptamers, i.e, those where the distance of Pt atoms 
belonging to two different clusters was less than two unit cells apart. We distinguish between 
two possibilities: 



1. Translationally Random Heptamers on Pt(lll) 

Since it is plausible that the orientation of the Heptamers is determined by the surface 
structure, we located each cluster Pt atom at the center of the triangles that make up the 
underlying hexagonal structure of the surface Pt atoms (i.e, a lattice site). This corresponds 
to a high barrier for rotation of a Heptamer-cluster as a whole. It leads to a 6-fold symmetry 
in the orientations of the Heptamers. Part of the surface is shown in figure ^ 



2. Translationally and Rotationally Random Heptamers on Pt(lll) 



Though the surface- induced orientation described in [II B 1| seems intuitively physically 
correct, we cannot exclude a completely randomly distributed orientation of the Heptamers, 
corresponding to a low barrier for rotation of a Heptamer-cluster as a whole. Therefore we 
investigated a second model, where also the orientation of the Heptamers is chosen randomly 
from a uniform distribution, as can be seen in figure 0. For both cases the coverage of the 
adsorbed Pt atoms was 5% of a monolayer. The He / Pt cluster+surface interaction for 
these models is given by eq.(|l]), with the exception of the second term which now contains 
an additional sum over all N Heptamers: 

N 7 

= E E VHe-pt{\v - r,,„|) + Vs{z) (5) 

n=l i=l 

Here rj „ = {xi^ni yi,ni Zi^n) represents the location of the ith Pt atom in the nth Heptamer. 
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C. Diffusion Limited Aggregates Grown by Kinetic Monte Carlo (KMC) 



We now proceed to consider not only compact clusters, but structures that can develop 
in kinetic conditions. The main factors determining the shapes of the resulting structures 
are: 

1. the deposition rate of the Pt atoms 

2. spatial migration of the Pt atoms upon the surface 

There is an extensive literature dealing with the modeling of the process of diffusion and 
aggregation upon a homogeneous surface^. The approach that will be taken here to sim- 
ulate the time-dependent structures obtained in the course of deposition and migration, is 
the KMC schemed. The essence of the approach is assigning probabilities for the hopping 
of adsorbed atoms between allowed sites upon the surface, which is assumed a lattice. To 
obtain probabilities for hops between the different configurations, assumptions about en- 
ergetics and interactions must be made. We shall present a kinetic model for growth of 
monolayer Pt islands on a flat Pt surface. Our purpose in examining this model was to 
provide realistic adsorbate structures for the He-scattering calculations. This enabled us to 
investigate possible ways of discriminating between these structures and formations of Hep- 
tamers. In order to produce realistic arrangements of adatoms on the surface we performed 
Kinetic Monte Carlo (KMC) simulations, using a model of pairwise additive interactions 
and a Time Dependent Monte Carlo (TDMC) procedure. The simulation took place on a 
hexagonal lattice of 100 ■ 100 unit cells with periodic boundary conditions. The simulation 
starts with two particles on the surface diffusing until a third particle adsorbs. This period 
is typically 0.05 sec (for simulating an adsorption rate of SOOsec/MzJl.) corresponding to 
about 10^ MC steps at room temperature. Then the three particles diffuse, a fourth one 
adsorbs, and so on until 500 particles, corresponding to 5% coverage, diffuse on the surface 
for 0.05 sec. The configurations attained are not at equilibrium. 
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a. The Interactions Model The simulations were performed using two different kinetic 
models for lateral interactions between the adsorbed particles. In the first model, every 
particle was assigned an energy according to the number of its nearest neighbors 



where Eq is the activation energy for diffusion at zero coverage, n is the number of nearest 
neighbors and e is the nearest neighbor interaction energy. Particle diffusion takes place via 
random walks between nearest neighbor sites. The hopping rate (transition probability per 
unit time) of a particle from site i, where it has rij nearest neighbors, to a neighboring site 
with Uf nearest neighbors, is taken as 



where ti^o = u exp[—EQ/ (kT)] is the hopping frequency of an "isolated" particle, tq = 
I/ujq, The average time interval between successive moves on the bare (zero coverage) surface, 
sets the unit of time for simulations performed at a given temperature T. Clearly, uji^f 
satisfies detailed balance, since Ui^f/uj^i = exp[—{Ef — Ei)/{kT)] with Ei = Eijii) and 
Ef = E{nf). Note that according to eq.(0), the transition rate depends only on the initial 
state. A possible dynamical interpretation of this model, originally suggested by Uebing and 
Gomerii, is that the rate of the transition from state i to state /, is governed by the rate 
at which the particle in state i escapes into a transition state (saddle point) which is the 
same for all possible transitions into different final states. Examples of structures produced 
by this model are shown in figure ^. 

In the second model we continued to use eq.(|^) for the transition rates, but now instead of 
eq.(^), the energy of a particle in a given configuration is: 



with n denoting, as before, the number of nearest neighbors of the particle under consid- 
eration, and s being the size of the cluster the particle belongs to. According to this model. 




(6) 




(8) 
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the nearest neighbor interaction energy e(s), depends on the cluster size (thus, indirectly, 
taking into account the effects of non-additive interactions). The assumption on the nature 
of the interaction is critical. It is clear that by assuming certain types of dependence of 
the microcluster energy upon shape and size, one can ultimately obtain any chosen type of 
kinetic behavior and global structure. We believe that the dependences we used, while def- 
initely not quantitatively valid, represent at least qualitative a plausible energetic behavior 
of a system of this type. The values of e(s) (for Pt atoms on Pt(lll)) were taken from the 
work of Rosenfeld et al.0, who used eq.(|D in their rate equation study of Pt diffusion and 
aggregation. More specifically, these authors used: 

e(s) = e(2) - (s - 2) Ae, with e(2) = O.SOOeV, and Ae = O.OldeV for 2 < s < 6, and: 
e(s) = e(6) - (s - 6) Ae, with e(6) = OMOeV, and Ae = 0.005eV for 7 < s < 14. Using this 
model in their rate equation approach, Rosenfeld et al. find a strong tendency for Heptamer 
formation at T > 4007^. On the other hand, our KMC simulations, using the same model 
for e(s), did not reveal a particular preference for Heptamer formation: mostly larger com- 
pact clusters were produced, with the number of constituent atoms varying from cluster to 
cluster. The difference must be attributed to the mean field type approximations inherent 
to the rate equation formalism. Having no reliable kinetic means at our disposal to produce 
compact heptameric islands, we used only the simpler first model to provide realistic ad- 
sorbate structures for our scattering calculations. The KMC simulations (for both models) 
were performed using the time-dependent-Monte-Carlo (TDMC) scheme, suggested inil. 
According to this scheme, instead of randomly chosing particles and accepting or rejecting 
moves according to the given transition probabilities (as in "traditional" MC simulations), 
one performs a move in any attempt, and propagates the time accordingly. More explicitly, 
we first calculate the average transition rate (r) = l/{uJi^f) out of state i, then randomly 
sample a given i —>■ f move with probability cjj^/ / J2 ^i^f s.nd, finally, perform this move 
and record the time elapsed as At = l/(r). Further details about this procedure can be 
found elsewhere@'ii. 
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D. Interaction Potential of He with DLA Structures 



We proceeded to perform scattering calculations from the DLA structures. Our purpose 
here was to check the degree of specificity of the features exhibited by the angular intensity 
distribution obtained from the Heptamer systems of sections [11 B 1| and 11 B 2 . The results 



of the KMC grown structures are the best model system we had at hand to simulate a 
"non-Heptameric" epitaxial surface. This model served as our control-case, against which 
the features identified for the Heptameric models were tested. We will consider more fully 
the question of what structural features can be extracted from the information contained in 
the scattering results from the DLA structures in a future article. The potential used for 
the He scattering calculations in this case was: 

N 

V{r) = Vs{z) + ^ VHe-Pt{\r - r,\) (9) 
1=1 

where Vs{z) and Vne-Pt are given by eqs. respectively, and ri now runs over the 

positions of all adsorbed Pt atoms. 

III. THE METHODS FOR THE SCATTERING CALCULATIONS 



We employed the sudden approximationlfjEj, which has proved very useful in studies 
of atom scattering from defects. Basically, the sudden approximation requires that the 
momentum transfer in parallel to the surface be small compared with the momentum transfer 
normal to the surface, i.e@ii, 

|K' - K| < 2h (10) 

where is the incident wave number in the z direction, K is the incident wave vector 
in parallel to the surface plane (zero in our simulations) , and K' is any intermediate or final 
wave vector in parallel to the surface plane which plays a significant role in the scattering 
process. For a surface area A of atom type X which contains a collection of defects of the 
same atom type, the angular intensity distribution is given by the function 
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2ir){p)+i p-(K'-K) 
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(11) 



where p = {x,y) and ri{p) is the scattering phase shift computed for fixed p, given in 
the WKB approximation by: 



oo 



Vip) = 1^^^^^ dz UK'- '^JIlX^lA -K\-K zt{p) (12) 



where ^t(p) is the classical turning point pertaining to the integrand in eq.(|T2]). Periodic 
boundary conditions were employed. Condition (|I0|) for the validity of the sudden approxi- 
mation is expected to break down for systems of high corrugation, and an isolated adsorbate 
on an otherwise flat surface generally represents, for realistic parameters, a very substantial 
local corrugation. Nevertheless, previous calculations have shown that the sudden approx- 
imation reproduces rather well many features of the scattering from isolated adsorbatesii. 
Features for which it breaks down are, e.g, intensity peaks due to double collision events (in 
which the incoming atom first hits the surface and then the adsorbate or vice versa), which 
are a particularly sensitive manifestation of a strong corrugation (we note that a double- 
collision version of the Sudden approximation has recently been developedH). The sudden 
approximation is thus useful for the interpretation of these features for which it works, and 
offers useful insights also at the points where it breaks down. The advantage offered by 
the simplicity of the sudden expression eq.(p!T|) is particularly manifest when considering a 
potential of the form: 

V{v) = V{z-i{^)) (13) 
This form implies that the classical turning surface is a translation of the shape function 

e(R)i: 

Zi(R) = e(R) + ;2o (14) 

Here measures the distance of the repulsive He-surface interaction from some reference 
plane. The choice of zq does not affect the physical observables (i.e., the intensities), zq is 



determined by the incidence energy and angle. Thus eq. ([l2|) can be rewritten as: 
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r/(R) = «-fc,(e(R) + ^o) (15) 
where a is a constant independent of R. Hence from eg. (pJ]) : 

thus avoiding the need even for the numerical evaluation of the phase-shift eq.(^^. It 
should be noted that within the sudden approximation the potential eq.(0) and a hard- wall 
potential with the same classical turning surface are essentially equivalent. Under certain as- 
sumptions concerning the interaction potential, the angular intensity distribution contains all 
the information about the scattering process. Indeed, it has been shown by Gerber et al.ilii 
that for an interaction potential of the form V{r) = Vq{z) + Vi{z) Q{x, y), if Vo{z) is known, 
it is possible to determine the form of Vi{z) and Q{x,y) from the angular intensity distri- 
bution. This was implemented with good accuracy for a problem where several diffraction 
peaks only are important, for a crystalline surface. However, in general it is a complicated 
experimental matter to measure the angular intensity distribution over the entire physical 
range, and one is usually restricted to several one-dimensional measurements. Interpreta- 
tion of the angular intensity distribution can very effectively be carried out by analyzing 
intensity maxima, which generally fall into three categories: (1) Rainbow peaks', related 
to scattering off inflexion points of the classical turning surface, (2) Fraunhofer peaks0, 
caused by interference of waves traveling around a spherical object, and (3) constructive 
interference (Interference) peaks, a consequence of interference of waves scattered from a 
periodic structure of some sort (which would be classified as Bragg interference for a fully 
periodic surface). The position of each peak, once identified, contains important structural 
information about the surface. Our ensuing analysis of the angular intensity distribution 
will proceed along these lines. Another, particularly simple feature from the experimental 
point of view, is the attenuation of the specular peak / by a defect. This attenuation is 
closely related to the cross-section of the defect. The cross-section may be extracted from 
eq.(|TTp via its experimental definitionil: 

/ = /oe-"^^ (17) 
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where / = Po^o, -^o represents the specular scattering intensity from the corrugated 
and flat (smooth, defect-free) surfaces respectively, n is the number of defects per unit 
surface area, and 9 their coverage. The parameter S in eq.(p!7|) can be interpreted as the 
cross-section for scattering by a single defect. When viewed as a function of the incidence 
energy, / oscillates with a frequency that is characteristic of the height of the defect above the 
surface. A Fourier analysis of these oscillations reveals the distribution of defect heights, and 
if performed at a wide enough energy range, will also reveal the z-direction "fine-structure" 
of these defects. This too will be an approach we shall implement in the next section. 
A combination of angular intensity distribution and / attenuation analysis proves to be a 
powerful tool in extracting information about defect structure. This will be developed in 
the course of our analysis of the results, later in this article. 



IV. RESULTS AND DISCUSSION 

The results will be outlined by considering separately each of the models described in 
section 0. Within each model, we will point out the main physical features present in the 
scattering calculations, pertinent to the classification of the type of disorder present on the 
surface. All scattering calculations were performed within the sudden approximation, with 
the He impinging on Pt(lll) at normal incidence. As will be demonstrated in detail in the 
next sections, there are very dramatic effects in the angular intensity distribution. For the 
case of the single adsorbed Heptamer, nearly all characterizing features of the Heptamer's 
geometrical structure can be discerned in an analysis of the intensity distribution. We are 
able to determine the symmetry, the linear extents, and the slope and number of inflexion 
points of the Heptameric shape function. Taken together, these define the Heptamer almost 
uniquely. For the case of a distribution of Heptamers adsorbed upon the surface, we see 
significant differences in the scattering results, between the model of translationally random 
Heptamers, and the model of rotationally random Heptamers. Both models, in turn, differ 
most considerably from the case of DLA structures. Furthermore, we are able to quantify 
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the differences among the various cases, relying upon the analysis of the single Heptamer. 
Thus we offer signatures of disorder structure based on the analysis of a He scattering 
experiment. We believe that the great majority of the pertinent features we observe will 
be experimentally observable. Much like in the case of scattering from ordered surfaces, 
an analysis of the angular intensity distribution, albeit considerably more complicated, may 
serve to almost fully classify the type of disorder present on a disordered surface. We proceed 
to demonstrate how such an analysis may be carried out. 

A. Physical Features in the Single Heptamer Intensity Distribution 

Figure ^ shows the results of calculations using the potential discussed in section p A . 
The collision wave-number of the impinging He was = 3.0bohr^^. The contour surface 
shown in figure ^ is the full 2D angular intensity distribution, i.e Pk^k' with the z-axis on 
a logarithmic scale. What is most striking about figure ^ is the unambiguous presence of 
structure in the intensity distribution: it has a sharp and well-defined angular appearance. 
The structure reflects the six-fold symmetry of the Heptamer, as indeed expected from a 
consideration of eq.(^). This property is easily shown to hold for other types of symmetry, 
i.e, the momentum-space view of the cluster via the angular intensity distribution retains 
the symmetry of the real-space cluster structure. This is a feature which should show up 
unmistakably in the appropriate experiment, if several on-plane directions are probed. Thus 
a small cluster on a surface has an experimental signature. 

1. Rainbow, Interference and Fraunhofer Effects 

In this section we will present in some detail an analysis of the major physical features 
present in the angular intensity distribution. This distribution is however much too detailed 
to give a full analysis. Such an analysis is definitely desirable to fully characterize the type 
of disorder present on the surface: in practice the angular intensity distribution contains 
enough information to leave no ambiguity. It will however be clear from the following, that 
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even a partial analysis may lead one a long way toward an understanding of the disorder 
one is dealing with. At the end of the section we summarize the main findings. 
Figure ^ shows a section along the axis (at ^Ky = 0) for scattering at a wave-number of 
kz = 3.0bohr~^, with a resolution of O.Obbohr"^ along the AK^ axis. The angular intensity 
distribution in figure | exhibits, in addition to a large specular peak, also a multitude of 
other peaks and "hidden" peaks, the pattern being obviously symmetric in the momentum 
transfer AK^ = K'^ — {K^ = 0) direction. In order to interpret the various peaks we 
proceed along the lines mentioned in the previous section, i.e we attempt to classify the 
peaks whenever possible as either Rainbows, Fraunhofers or Interferences. 

a. Rainbows: The take-home message from the ensuing analysis, is that Rainbow peaks 
appear as a dominant and experimentally observable feature in the angular intensity dis- 
tribution of He scattered from a single compact cluster, and that their position contains 
critical information about the surface morphology in the vicinity of the cluster, and about 
the interaction potential with the He. Measurement of the Rainbow peaks alone, is sufficient 
to calibrate important points in a simple model potential, and is hence a powerful tool in 
analysis of the surface disorder. Following will be a brief reminder of the origin and physical 
significance of Rainbows (this subject has been treated in much detail elsewhere; see, e.g.,ii), 
and a detailed analysis of Rainbows-related features in our scattering results. 
It is useful to employ a stationary phase, approximate evaluation of the sudden scattering 
amplitude, as i n@. The crude quasi-classical stationary phase hmit of eq.(|ll|) predicts a 
Rainbow whenever the phase-shift has an inflexion pointil: 

ri"{x) = (18) 

The Xq satisfying this condition dominates the scattering. In the classical limit this 
should show up as a singularity in the angular intensity distribution at momentum transfer 
AKx such that 

AK^ = -2r]'{xo) (19) 
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where xo is the stationary phase point of the integrand eq.([TT|) for momentum transfer 
AKr^. The singularity of the crude classical evaluation is smoothed into a finite peak in the 
more refined sudden calculation. Figure ^ shows the phase-shift function r]{x) at ?/ = 0, along 
with its inflexion points for kz = 3.0 bohr~^. Due to the inversion symmetry of the Heptamer 
about the y-z plane, the inflexion points appear as symmetric pairs along the x-axis. Each 
such pair gives rise to corresponding symmetrically located Rainbow peaks, as is clear from 
eg . ([19D . The similarity between rj^x) and the Heptamer's isopotential (flg. [T^ ) is understood 
from eq. ([l^) . Based on this similarity, we can identify the physical origin of the Rainbow 
peaks: while it may be difficult to demonstrate this for a general potential (but immediate 
for, e.g., eq.(|13])), it is clear that in the present case the inflexion points of the phase shift 
correspond uniquely to those of the turning surface. From flgure ^it can be seen that at most 
four distinct Rainbows peaks are expected, one for each inflexion point. The points numbered 
1 and 2 are due to the Heptamer's top corrugation, number 3 due to its side, and number 
4 due to the corrugation induced by the long-range attractive part of the He-Pt interaction 
at the vicinity of the surface. We expect such inflexion points (and the resulting Rainbows) 
to be widely present for adsorbates on flat surfaces, regardless of the quantitative details of 
the potentials. In our case, the corrugation associated with point 4 is too small to enable 
R4 to be discerned from the specular peak. R3 is clearly associated with a double collision 
in which the He flrst strikes the Heptamer's side and then undergoes a second collision 
with the flat Pt surface. As mentioned before, such Double- Collision Rainbows cannot be 
described within the sudden approximation, because it does not include double-collision 
effects. However, using the quantitative predictions of eqs.(|18|,|19D, we identify the peaks 



at AKx = ±0.33 bohr^^ and the hidden peaks at AKx = ±0.13 bohr^^ as Rl and R2 

respectively. The peaks have an intensity of 10~^ and 5 ■ 10"'^ respectively with respect 
to the specular, and should be well within experimental reach (the intensity predictions of 
the sudden approximation for RB peaks are generally in good agreement with the exact 
quantum-mechanical calculations; see e.g.ii). We did not flnd a point Xq that satisfles the 



Rainbow condition eq.(^), i.e, inflexion points of the phase shift, for the maxima at other 
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/S.Kx values. We found it instructive to compare our results with those obtained by scattering 
from the same system, at a significantly lower incidence energy. Our motivation was the 
expectation that at a lower energy, the He probes the Pt system at a more global scale, and 
is hence more sensitive to the overall structure of the defect and less to local details. This 
follows also from the reduced penetration depth, resulting in a reduced effective corrugation. 
The scattering calculations were repeated for the He impinging on the Pt surface at normal 
incidence, with a wavenumber of kz = 2.0 bohr~^. Figure |^ shows the same section as figure |^ 
for this energy, with a more than 4-fold improved resolution of 0.012bohr~^ (compared to 
the scattering at kz = S.Obohr"^). Figure ^ shows the corresponding phase-shift at y = 0, 
along with its inflexion points. The inflexion points at the Heptamer's top do not persist, 
as at this low incidence energy, the He particle sees a higher Coulomb barrier in between 
the top Pt nuclei. Hence it is less effective in penetrating this "soft" region, resulting in 
the smoothed appearance of the Heptamer's top. The absence of inflexion points does not, 
however, rule out the existence of Rainbows: it must be remembered that conditions (|T8D 
and (|r^) follow from a crude stationary-phase analysis and are to be considered merely as 
sufficient conditions for Rainbows. We do however expect a significant reduction in the 
intensity of the Rainbows. Furthermore, analysis of scattering results at kz = 4.0 and 
5.0 bohr'^^ (not presented here) reveals that the position of the Rainbows is nearly constant, 
showing only a slight distancing from the specular peak with increasing kz, as expected from 
the corresponding increase in corrugation. Returning to figure |^ we thus identify the peak 
at AKx = 0.15bohr^^ as R2. Compared with figure |^, its intensity is seen to have indeed 
decreased by as much as two orders of magnitude. Due to the improved AK^ resolution we 
are now also able to discern the hidden peak at AK^ = O.OlSbohr^^ as R4, resulting from 
the long-range attractive part near the edge of the Heptamer. Rl is missing, apparently 
masked by other peaks. 

b. Interference Peaks: The origin of these peaks is the constructive interference between 
He waves scattered from different atoms. On a perfect crystalline surface, these would be 
referred to as Bragg interference peaks. Here, we expect them to arise from constructive 
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contributions from neighboring Pt adatoms in the Heptameric cluster. No contribution 
can be expected from the surface, as it was treated as completely flat. This is a realistic 
assumption for Pt (111), which has an extremely smooth surface, whereas the attraction of 
the Heptamer's top electrons to the underlying surface will result in a significant corrugation. 
Interference peaks are expected at Aii'^. = a = 5.23 bohr being the Pt unit cell distance, 
n an integer, for scattering along the x-direction. However, it turns out that even though 
these peaks exist, they are of negligible magnitude (< 10~^ with respect to the specular) 
and not within experimental realm. We shall return to this point when we discuss the more 
realistic system of the next section. 

c. Fraunhofer Peaks: Finally we attempted to assign the remaining peaks using the 
Fraunhofer diffraction model for scattering from a hard hemisphere on a fiat, hard surface, 
given by Lahee et al.0. The Fraunhofer diffraction intensities given by this model 



are: 



P{e) oc \ Ji{dk sm{e)) csc(e)r (1 + cos{e)Y (20) 

where 6 is the sum of the incidence and scattering angles (simply the scattering angle 
in our case), d the hard-sphere radius, k the collision wave-number, and Ji denotes the 
Bessel function of 1'^* order. Given that eq.(|20) is derived under assumptions valid for the 



scattering of light waves on small obstaclesE3, we did not expect it to accurately predict the 
intensities of the Fraunhofer peaks (see also ref. 0). Hence rather than trying to fit eg . (^OD 
or its parameters to our data, we used a simple prediction resulting from eq.(pOD: that the 
spacing between Fraunhofer peaks as a function of = k sin (6') is very nearly constant. We 
expect this property of the Fraunhofer model to be more robust than the intensities, as the 
latter are modified by modulations due to other sources (such as Rainbows). While constant 
spacing cannot be claimed to be unique to the Fraunhofer model, we are not aware of any 
other relevant model with the same property, which is in addition appealingly simple. Thus 
we tried to fit a linear regression for the spacings of the remaining peaks. This could be done 
with a maximum error of 6% per peak and a coefficient of determination = 0.995 (fig. 
H). The constant peak-spacing prediction is expected to persist along the principal axes 
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also, with a more sophisticated two-parameter model of scattering from a hard ellipsoid 
replacing the hemisphere. Indeed, we find a value of AKy{F) / AKx{F) = 1.16 for the ratio 
of the predicted Fraunhofer peak spacings in the AKy and AK^ directions, compared to 
a geometric ratio of R^/Ry = 1/ cos(7r/6) = 1.15 between the major and minor axes of 
the Heptamer (the inversion is expected due to the conjugate relationship between real and 
Fourier space). We thus confidently identify the remaining peaks as Fraunhofer diffraction 
peaks. Once again, it was instructive to compare these results to the case of a lower incidence 
energy. Most of the structure in |^, at kz = 2.0bohr~^, is in fact due to Fraunhofer peaks. 
Indeed, the very regular appearance of the non-rainbow peaks and their almost smooth 
envelope is very reminiscent of a Fraunhofer interference pattern. Relying on the constant 



spacing property of eq. (pOD again, we were able to fit a linear regression to the peak spacings 
with a maximum error of 1.3% per peak and a coefficient of determination R^ = 0.99995 
(fig. |]). The dominance and high regularity of the Fraunhofer pattern at this energy is 
explained by the degree of smoothness of the Heptamer itself, in addition to it being a single 
corrugation on an otherwise flat surface. As we shall see in the next section, this feature 
breaks down when many adsorbates are present on the surface. The Fraunhofer pattern 
also obscures any Interferences, which were in any case expected to be small due to the low 
corrugation of the Heptamer's top. 



2. The Cross-Section of a Single Heptamer 

The cross-section, being the simplest parameter to extract experimentally (indeed, the 
bulk of the experiments done so far in He scattering have been cross-section measurements. 
See, e.g, ref.il), is probably also the single most descriptive parameter of the surface disorder. 
As we shall see in this section, the cross-section has an oscillatory behavior in the incidence 
energy, and there is a mathematical procedure, outlined in appendix A, which enables a 
determination of the cluster heights distribution from these oscillations alone. However, it 
should be realized that the necessary large scale variations of the incidence energy are not 



22 



easily realizable experimentally. 

Figure |10| shows the cross-section for He scattering at normal incidence of Pt Heptamers 
at energies corresponding to kz = 1, . . . , 4 feo/ir^^. The cross-section was calculated from 
eg. (p!7D , which reads: 

.^WoZ (21) 



n9 

when inverted for S. In our case n6 = 1/ {total area) (there is one defect - the Heptamer 
- on the entire surface) and Jq = 1 (the defect-free surface is flat). Two main points are 
worth noting from figure |TU|: 

1. The value of the cross-section, which lies between 320 and 430 A^: 

The cross-section of single Pt adatoms on Pt(lll) has not been measured yet. We 
assume it to be similar to that measured for a number of adsorbed molecules and for 
vacancies, i.e about 120 for a room temperature He beam incident at Oi^c = 40°ii. 
This then yields an effective Pt radius of 9A, and, using the superposition principleii 
together with a Pt unit cell distance of 2.77 A, we arrive at a cross-section of about 
7r(2.77-|-9)^ cos(40°) = 195A^ for He scattering off a Pt Heptamer at normal incidence. 
Thus our calculated cross-section is about twice as large as the empirical estimate. This 
almost certainly reflects the inaccuracy in our value for cg, as explained in section [11 A. 



However, the actual cross-section of a Pt atom on Pt(lll) still remains to be measured 
and it may come out higher than the value employed above. At any rate, our main 
purpose in this article is to provide qualitative insight into the scattering features. We 
thus note that it is important to obtain a reliable experimental value for the single 
adatom cross-section, which once known, can be used to estimate the cluster cross- 
section. This in turn, can narrow down strongly the number of possible geometric 
configurations to be considered in attempting to identify the cluster type. 

2. The cross-section oscillations: 

These arise from interference between waves scattered from the flat part of the Pt 
surface and the Heptamer's top. Clearly the period of oscillations is related to the 
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height difference between the surface and the top of the corrugation. A Fourier analysis 



of the oscillations in figure |T0| is shown in figure |TT|, and yields a height difference of 
about 5.5bohr. This compares quite well with a direct calculation from the classical 
turning surface (see figure |T2|), which yields an average height of about 6.5bohr. In 
the Heptamer case there are effectively two heights to be considered (marked hi, /12 in 
figure [12|) . To see the "fine structure" (/12) from the cross-section oscillations, clearly a 
large range of incidence energies must be employed. In the case of the Heptamer, our 
calculations show an unrealistic several hundreds of eV would be required to discern 
h2. However, with adsorbates composed of several atom types, or even an adsorbate 
of different substance than the surface (hke in the adsorbate-induced restructuring of 
surfaces of Cu(110)-(2 x 1) reported by Kern et al.0), the situation could be quite 
different and several experimentally observable heights could exist within one cluster. 
A theory explaining how to obtain the distribution of cluster heights from the cross- 
section oscillations is hence desirable, and appendix A provides a detailed calculation 
of this type in the case of aV = V{z — ^(R)) potential. 



3. Summary of Results for Single Heptamer 
We sum up now briefly the main findings of this section: 

1. Rainbow, Fraunhofer and Interference peaks are all expected to be seen with varying 
intensities in scattering from a Heptamer in the same collision conditions (i.e the same 
experiment). 

2. The Rainbow peaks due to the Heptamer's top corrugation are expected to be the 
most pronounced feature and to lie closest to the specular. Fraunhofer peaks, probing 
the size of the adsorbate, will vary in intensity and the first few may be masked by the 
much stronger Rainbows, however they are expected to be found at nearly constant 
spacing. Interference peaks due to constructive interference between waves scattered 
from the Heptamer's top atoms may be too weak to be measured experimentally. 
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3. As the incident He's energy is lowered, we expect the angular intensity distribution 
to become more regular, with the Praunhofer peaks gradually overtaking the Rainbow 
peaks in dominance. This reflects the prevailing of the small-scale corrugation at high 
energies, which is smoothed out when the De Broglie wavelength of the incident He 
becomes more comparable with the Heptamer's dimensions at < 1.2bohr~^. 

4. The cross-section is a most useful tool in that it is both very simple to obtain exper- 
imentally and also provides very important information about the defect's size and 
topography. The size information can be used to narrow down the possible candidates 
for cluster shapes, and the topography provides an excellent characterization of the 
defect's appearance. 

5. As each of the scattering features discussed here is sensitive to different aspects of the 
interaction potential, experiments measuring all these features should provide a nearly 
complete signature of a Heptamer, if present on a surface. This conclusion should hold 
quite generally for other adsorbate shapes. For example in the Heptamer case, the 
prevailing of Fraunhofer peaks at low incidence energy suggests one is dealing with a 
nearly hemispherical, i.e. compact object. One may then combine the measured cross- 
section and the diameter estimated from the Fraunhofer peaks' position to calculate 
how many Pt atoms participate in a given cluster. The number of observable Rainbow 
peaks gives one a lower bound on the number of distinct inflexion points in the cluster, 
which can be used to geometrically reconstruct the cluster. The positions of the 
Rainbow peaks tell one about the degree of corrugation of the cluster. An analysis 
of the cross-section oscillations reveals the number of different heights in the cluster, 
their values and relative frequency of appearance. Taken together, these data combine 
to determine the cluster quite fully. 

We proceed now to examine the more realistic model of scattering from a collection of 
randomly adsorbed Heptamers. 
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B. Scattering Results from Randomly Adsorbed Heptamers and DLA Structures 



In the following sections we analyze in some detail the results for scattering from the 
randomly adsorbed Heptamers and the DLA structures. The main message here, is that 
structural features persist in the angular intensity pattern of distributions of compact clus- 
ters, and can be quantified to extract interesting parameters of the disorder. Very little 
structure, if any, however, persists in the intensity distribution of DLA, in sharp contrast 
to the compact cluster case. Thus a scattering experiment can readily distinguish between 
compact clusters and DLA adsorbed on a surface. 

We proceed to analyze our results. A summary is presented at the end of the section. 



1. Translationally Random Heptamers on Pt(lll) 



Figure 13 shows the results of the scattering calculations from the system described in 



section |II B 1|, again for kz = 2.0bohr ^. There is a rich structure in this angular intensity 



distribution and we shall not attempt to assign every peak. However, relying on the analysis 
of the previous sections, we identify Rainbow peaks at AK^ = ±0.024 6o/ir~^ (1^4), AK^ = 
±0.12bohr^^ {R.2), and AK^ = ±0.36 6o/ir~^ (R-l)- These values agree very well with 
those found in the single Heptamer case, and we attribute the small discrepancies chiefly 
to numerical deviations. The reversal in relative intensities of R4 and R2 (R2 appears 
more intense now) is believed to be the result of the different role played by the Fraunhofer 
diffraction pattern in the current system. The optimal conditions for Fraunhofer diffraction 
are when there is no interference between neighboring clusters, as was certainly the case for 
the single Heptamer. In the present system the average distance between Heptamer centers 
is rather small, giving rise to significant interference between neighboring Heptamers. As 
a consequence the Fraunhofer pattern is strongly attenuated. The increased intensity of 
R4 in figure ^ is hence attributable to an underlying Fraunhofer peak. Considering the 
shallowness of the attractive well giving rise to R4 (see figures |],fig:PS2), the "normal" 
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state of affairs is restored in figure |T3|, where R2 is the more intense. The attenuation 
of the Fraunhofer pattern allows for the reappearance of Rl, which was missing (masked) 
in figure |^. The strong Interference peaks at Ai^x = ±1.2bohr~^ can of course not go 
unnoticed. Their intensification must be attributed to the constructive contributions of all 
Heptamers, combined with the diminished masking by the Fraunhofer peaks. 

2. Translationally and Rotationally Random Heptamers on Pt(lll) 



Figure 14 shows the results of the scattering calculations from the system described 



in section PB 2|, at the same incidence wavenumber of = 2.0bohr ^ as in the previous 



section. It is most instructive to compare figure |1^ to figure |13[ Clearly the envelope of 
the intensity spectrum is much reduced, as expected from the disappearance of most of 
the constructive contributions when the Heptamers assume random orientations. Yet the 
Rainbow peaks persist, albeit at a lower intensity, and we again identify R4 at AK^ = 
±0.02Abohr~\ R2 at AK^ = ±0.12bohr-^ and Rl at AK^ = ±OMbohr'^. Like the 
Rainbows the Interferences at AKx = ±1.2bohr^^ are also much weaker, corresponding to 
the fact that there are now fewer constructive contributions in the original AK^ direction, 
and most of the intensity is present in other directions. 

3. Scattering Pattern for DLA Structures 

The scattering results from this system, described in section |II C| (first model), are pre- 
sented in figure |15[ The incident He's wavenumber was once again = 2.0 bohr"^. The most 



prominent feature in figure [T^ as compared to figures ^ and |14| is the almost complete dis- 
appearance of structure in the angular intensity distribution. The peaks that can still be dis- 
cerned are the Interferences at AK^ = ±1.2bohr~^ and a new set at AK^ = ±0.096bohr~^ . 
The existence of the small Interference peaks is explained by the hexagonal unit cell struc- 
ture of the islands. This is of course far from a complete hexagonal surface at the 5% 



coverage of figure |T^, but is sufficient to produce some constructive interference between 
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neighboring Pt atoms. The peaks at = ±0.0966o/ir~^ are not at the position of any of 
the Rainbow peaks from the previous two sections, but can still be considered as Rainbows, 
this time an Rl or R2 (they are identical from symmetry considerations except for at the rim 
of the islands). They arise from the weak corrugation at the top of the monolayer islands. 
Another dominant feature of figure is the broad specular peak. This means that intensity 
is distributed almost continuously into nearly specular directions. We shall not attempt an 
analysis of the width of the specular peak and its behavior as a function of AK^ here. We 
shall limit ourselves to mentioning that these features are related to the fractal dimensions 
of the Pt islands^. A detailed examination of this issue in the context of He scattering 
will be presented in a future article (see also ref. 0). Needless to say, no Fraunhofer peaks 



pattern is observable in figure 15. 



4- Summary of Results for Randomly Adsorbed Heptamers and DLA Structures 
To summarize our findings in the last sections: 

1. There is a rich structure of several peaks in the case of He scattering off randomly 
adsorbed Heptamers. Much of this structure persists, albeit at a lower intensity, for 
the case of randomly rotated Heptamers. However, in the case of randomly structured 
islands, there is nearly no trace of structure in the angular intensity distribution. 
Thus a powerful preliminary criterion to distinguish structured clusters from randomly 
shaped ones, is the respective structure of their angular intensity distribution. 

2. The source of the structure in the angular intensity distribution, in the compact- 
cluster case, are primarily the Rainbow peaks. These are at the same time very robust 
to perturbations of the structure (e.g a random rotation of all clusters) and strongly 
position-related to the fine details of the cluster structure. Thus the Rainbow peaks 
may serve as a crucial component in the identification of surface structure. 
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3. At this relatively low concentration, the Praunhofer peaks are of secondary importance. 
They are quite sensitive to adsorbate proximity and shape. 

4. Interference peaks tend to reflect the underlying surface symmetry as displayed by the 
adsorbed clusters and are hence not of major help in determining adsorbate structure. 

5. A broad specular peak without much other structure is indicative of randomly shaped 
islands on the surface. 

6. The question arises whether the identification of cluster type is still as clear as in the 
single Heptamer case. We believe the answer to be affirmative. The major differ- 
ence between the single Heptamer and randomly adsorbed Heptamers cases was the 
prevailing of Fraunhofer peaks in the former. This resulted from the lack of inter- 
ference between neighboring clusters in the single Heptamer case. This situation can 
however be experimentally approached by reducing adsorbate concentration. If com- 
pact clusters are formed, they will then clearly create a Fraunhofer diffraction pattern. 
With increasing concentration interference between neighboring clusters will destroy 
the Praunhofer pattern, but if the average cross-section remains the same, one may 
conclude quite safely that compact clusters remain present. If compact clusters do 
not form at concentrations where there is negligible interference between neighbor- 
ing clusters, one must resort to relying on Rainbows and cross-section analysis alone. 
However, as discussed in the single Heptamer case, this analysis can lead one a long 
way toward a cluster signature. 

V. CONCLUDING REMARKS 

In this study we focused on the utility of TEAS as a tool for learning about surface 
disorder. More specifically, we sought to demonstrate that the angular intensity distribu- 
tion of He atoms scattered from adsorbed clusters contains information which can be used 
to determine most, if not all the characteristics necessary for a signature of a cluster. This 
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determination is based on the identification of Rainbow and Praunhofer peaks in the angular 
intensity distribution, and an analysis of cross-section oscillations dependence on incidence 
energy. All these features, by probing different aspects of the interaction of the adsorbed 
clusters and the colliding atoms, when taken together, combine to provide an almost com- 
plete identification of the cluster. Although our results were obtained for the particular 
case of Heptameric clusters, common physical sense dictates that the qualitative features 
should persist for other cluster types. That is, as we have demonstrated for Heptameric 
clusters, sharp and well-defined characteristics should be extractable from the angular in- 
tensity distribution. This should certainly be so for the case of a single cluster (or a dilute 
concentration of clusters), and with some more labor for the case of a distribution of clus- 
ters. For example, we expect that in the case of a distribution of compact trimers, only 
one Rainbow peak will be seen, the Fraunhofer peak spacing will reflect the structure of the 
underlying triangular base, and the cross-section will indicate the presence of three atoms. 
Thus the scattering can serve as a crystallographic tool for disordered surfaces. The major 
reservation regarding the last statement is, that a distribution of sizes and shapes of clus- 
ters, may obscure the otherwise sharp features present in the diffraction spectrum. Viewed 
from another perspective, as we have shown for the case of translationally random Hep- 
tamers versus translationally and rotationally random Heptamers, versus DLA, inspection 
of the angular intensity distribution is a means to distinguish one type of disorder from 
another. We are aware that the effects discussed here are hardly easily detectable. Yet, we 
have performed our calculations in a realistic energy regime, and by our estimates most of 
them should be observable with present day techniques and resolutions. On this basis one 
may anticipate considerable advance in knowledge of compact clusters adsorbed on surfaces 
and their interaction with incoming atoms, once the pertinent experiments are performed. 
An issue we have however only scratched the surface of, is that of scattering from DLA. 
Whereas we have shown that the structure of the angular intensity distribution in this case 
is substantially different from the compact clusters case, and hence such structures are easily 
distinguished from compact clusters, it is as of yet an open question what information can be 
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extracted from the angular intensity distribution in this case and how this information could 
be used to characterize such structures. Work on this issue is in progress and once completed 
should significantly enhance our ability to understand surface clusters. In light of the very 
impressive results obtained recently by means of the STM technique (see, e.g.,0S), a few 
comments on the relative position of the two techniques are in order. In our view, TEAS and 
STM are really complementary techniques, in the sense that STM probes the local structure 
(several millions of A^), whereas TEAS probes the global structure. In addition, TEAS is 
better at probing defects with high mobility, and does not damage the surface. However, 
there is another significant difference, which is frequently overlooked. This is the fact that 
STM mostly measures the electronic density near the Fermi level, whereas He scattering 
is sensitive to the total electronic density, via the He-surface interaction potential^. Thus 
although He scattering suffers from the significant disadvantage of displaying the surface in 
momentum space, whereas STM shows the surface in real space. He scattering is applicable 



applied only to metal surfaces. We hope to have convincingly demonstrated that even in a 
field in which STM has taken a definite lead. He scattering will have much to offer once the 
pertinent experiments are performed. 
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Appendix 

Calculation of Surface Heights Distribution from Specular Peak Attenuation: 



It suffices to measure the intensity of the specular peak over a range of incidence energies to 
obtain the distribution of surface heights. Consider eg. (piBD , which was obtained under the 
assumption of a potential of the form (eq.(|l3D): 



V{r) = V{z-aR)) 



Assuming a specular collision (K' = K), we expand eq.(|T6]) to obtain: 

loih) = Pk^k = -^jj dRdR'e^"'^^^^^''>'^^^'>^ (22) 

We know that the cross-section, and hence the specular peak intensity, oscillates with 
the incidence energy. To see this, we Fourier transform lQ{kz). Defining first the difference 
between shape functions: 

e(R,R0 = -2[e(R) -e(RO] (23) 

the Fourier transform is given by: 

I{z) =J^[Io{h)] = 1. I I dRdR' I (ifc^e'^^^t^-^^^'^')! 

= ^ J J dRdR' S[e{R,R') - z] (24) 

Equation (p^ ) can in principle be solved for T{z), given the shape function ^(R). We 
claim that its normalized value for a given height above the surface z, is the relative frequency 
of appearance of that height. In other words, I{z) is the distribution of surface heights. To 
see this, consider the following example, corresponding to a radically simplified "lego" -model 
of 3D epitaxial growth. Suppose that atoms can adsorb only along the x-axis, and only one 
nucleation-site is available. Suppose further that the resulting defect's shape-function is 
adequately represented by a stacking of rectangles of width w and height h. The situation 

35 



is represented schematically in figure ^ Let the number of atoms making up the base of 
the defect, i.e, the total number of heights be N. We denote the height at Xj by hj. Let the 
number of distinct heights be M, i.e, there are M different heights {h}fLi^ each appearing 
m times. Thus: 



r M 



w 



Next define a height-function: 



A,-(x) 



ii X < jw OT X > {i + l)w 
—2hj ii j w < X < {j + 1) w 

Using this, we have for given x': 

Aj{x) - Ay{x') i{0<x,x' <L 

if < X < L and {x' < ot x' > L) 



Q{x — x') = < 



A,(x) 







Now we can evaluate Ell 



if < x' < L and (x < or x > L) 



if X, x' < or x,x' > L 



dx I dx 5[Q{x,x') — z] = dx' [h + I2 + I3 

J J —00 



where: 



Jo 



dx 6 



Aj{x) - Ajv(x') ~z if < x' < L 



Aj{x)-z 



if x' < or x' > L 



(25) 



(26) 



(27) 



(28) 



(29) 



dx 6 



-Aj>{x') -z if < x' < L 



if x' < or x' > L 



(30) 



dx 6 



-Af{x') -z if < x' < L 

—z if x' < or x' > L 



Expansion of each term yields: 



(31) 
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Ill 



j £ dx' £ dx 5[ Aj(x) - Ajv(x') - z\ + 



II 



dx' r dx5[-^j\x') - 

Jo J-oo 



z] + 



dx' dx5[-/^j,{x')- z] 



II 



dx' dx5[/^j{x) - z] + / dx' dx5[-z] + 

J J—oo Jo J —oo J —oo 



/O /■oo 
dx' / dx S[ —A 
-oo J L 



/roc rly roo rO 

dx' J dx5[Aj{x) - z] + dx' J dx5[-z] + 



roo roo 

J dx' J dxS[—z] 

It is simple to check that all terms overbraced with the same numeral are indeed 
Their evaluation yields: 

l = {M-^)\{z) 
N-i ,r-.-4-ii,„ T N-i 



roo J' / r(j+i-)w r 

II = / dx'Y. dxS{-2 hj -z) = {M--)wJ2 '^(2 hj + z) 



N-l N-1 rij+l)w 

^ ^ / dx' dx d{-2 hj + 2 hj, - z) 



N-l N-l rij+l)w 
III = ^ ^ 

j'=0 j=0 

N-l N 

j'=0 j=0 



The defect-free zone is of length: 



Lo^2(M-^) 
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Collecting all terms and inserting the definitions p5| ) and (|38|), we obtain finally the 
normalized result: 

1 2 1^^ 

I{z) = 5{z) + ^ (^) E ^(^ + 2 /.) + 

I' „ ■' 

IV 

1 T 2 M M 

W^T^ EE^^'^^'^[^-2(^.'-/.)] (39) 
III' 

The physical significance of the terms in eq. ( |39|) is as follows: 

1. I' represents the contribution of the flat part of the surface: a unit-sized peak at ^ = 0, 
the reference plane. 

2. II' is the most interesting term: it contains a contribution from all heights present on 
the surface k, each weighted by its frequency of appearance rij. The sum over these 
contributions is the sought after distribution of surface heights. It is scaled as 

3. Ill' is a sum of peaks at all possible height differences. This information is redundant 
as it can be extracted from II'. It is scaled as (-^)^ (-^)^ and will therefore be negligible 
compared with II'. 

Equation (^) is a most useful result in our opinion. It shows that a very significant 
statistical feature of the disorder, the distribution of surface heights, can be extracted with 
ease from a He scattering experiment. Experimentally, all that needs to be done is to 
measure the specular scattering at a wide incidence energy range. A Fourier transform will 
reveal the desired distribution. 

The situation for the more realistic situation of defects on a surface (as opposed to a line) 
is even better: our calculations show that the "noise" term III' is scaled down by a factor 
of (-^)^ with respect to the "signal" term II'. 



38 



FIGURES 

FIG. 1. Schematic representation of translationally random Heptamers at 5% coverage on 
Pt(lll). Although the Heptamers occupy random lattice sites, they clearly have orientations 
compatible with the underlying hexagonal symmetry of the surface. 

FIG. 2. Translationally and rotationally random Heptamers at 5% coverage on Pt(lll). Here 
the Heptamers do not obey an orientational constraint. Instead their orientation is picked at 
random from a uniform distribution. Their centers of mass occupy random lattice sites. 

FIG. 3. Ft islands produced by KMC simulations, at 5% coverage on surfaces of 100*100 unit 
cells, according to the model with energy independent of cluster size (first model; see text). The 
surface is a parallelogram with an angle of 60°. Periodic boundary conditions are employed, so 
straight edges usually imply that the cluster is continued at the opposite boundary. The configu- 
rations displayed differ only in the random number sequences used to generate them. 

FIG. 4. The full angular intensity distribution for He scattering from a Ft Heptamer on a Ft 
surface, at kz = 3.0bohr~^ . Axes: AK^ (horizontal), Ai^y (vertical), log intensity (out of page). 
The six-fold symmetry of the Heptamer is clearly maintained in momentum space. 

FIG. 5. Log He scattering intensity as a function of momentum transfer in x direction, at 
kz = S.Obohr^^ . Positions of rainbow peaks are indicated. 



FIG. 6. Phase-shift function r]{x) as defined in eq.(12), at kz = S.Obohr ^: section along x 
axis. The inflexion points are indicated. 

FIG. 7. Log He scattering intensity as a function of momentum transfer in x direction, 
at kz = 2.0bohr^^ . Positions of rainbow peaks are indicated. A strong Fraunhofer interference 
pattern is visible. 
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FIG. 8. Constant spacing between peaks in Praunhofer model: the positions of peaks identified 
as Fraunhofers in the scattering calculations are displayed as a function of peak number, for two 
incidence energies. The (nearly) constant spacing between peaks is predicted by eg. (pO|) . 

FIG. 9. Phase-shift function r]{x) at kz = 2.0bohr^^: section along x axis. The inflexion points 
are indicated. 

FIG. 10. Total cross-section (in A^) of the adsorbed Heptamer, as a function of He incidence 
wavenumber kz{bohr~^). 

FIG. 11. Fourier transform of the specular intensity versus the height above the surface. The 
peak at z ~ 10.56o/ir corresponds to a distance of ~ 5.5bohr from the Heptamer's top to the 
surface, as explained in appendix A. 

FIG. 12. Heptamer classical turning surface at y=0 for two He incidence wavenumbers: 
kz = 4:.0bohr~^ and k^ = 1.7bohr~^ . The two heights relevant to the distribution of surface heights 
are indicated. 

FIG. 13. He scattering intensity distribution from translationally random Heptamers, along 
AKx[bohr^^] axis. He incident at kz = 2.0bohr~^ . Rainbow (R) and constructive interference 
(Int.) peaks are indicated. 

FIG. 14. He scattering intensity distribution from translationally and rotationally random 
Heptamers, along AKx[bohr^^] axis. He incident at kz = 2.0bohr~^ . 

FIG. 15. He scattering intensity distribution from DLA, along AKx[bohr^^] axis. He incident 
at kz = 2.0bohr-^. 

FIG. 16. "Lego" model of 3D epitaxial growth. Schematically drawn is the shape function of 
the surface with a single defect adsorbed upon it. 
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